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Abstract 

We consider a processor sharing storage allocation model, which has 
m primary holding spaces and infinitely many secondary ones, and 
a single processor servicing the stored items (customers). All of the 
spaces are numbered and ordered. An arriving customer takes the 
lowest available space. We define the traffic intensity p to be A//i 
where A is the customers' arrival rate and \i is the service rate of the 
processor. We study the joint probability distribution of the numbers 
of occupied primary and secondary spaces. For < p < 1, we obtain 
the exact solutions for m = 1 and m = 2. For arbitrary m we study 
the problem in the asymptotic limit p f 1 with m fixed. We also 
develop a semi-numerical semi-analytic method for computing the 
joint distribution. 
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1 Introduction 



We consider the following storage allocation model. Suppose that near a 
restaurant there are m primary parking spaces and across the street there 
are infinitely many additional ones. However, the restaurant has only one 
waiter who serves all of the customers. All of the parking spaces are num- 
bered and ordered; the one with rank = 1 is closest to the restaurant and 
the primary spaces are numbered {1, 2, 3, , m}. We assume the follow- 
ing: (1) customers arrive according to a Poisson process with rate A, (2) the 
waiter works at rate /i, (3) an arriving car parks in the lowest-numbered 
available space, and (4) if there are N customers in the restaurant, the 
waiter serves each customer at the rate n/N. This corresponds to a pro- 
cessor sharing (PS) service discipline. 

Dynamic storage allocation and the fragmentation of computer mem- 
ory are among the many applications of this model. We define N\ to be 
the number of occupied primary spaces and N 2 to be the number of occu- 
pied secondary spaces. Then we define S to be the set of the indices of the 
occupied spaces, and the "wasted spaces" W are defined as the difference 
between the largest index of the occupied spaces (Max S) and the total 
number of occupied spaces (|S| = Ni + N 2 ). Coffman, Flatto, and Leighton 
showed that for the processor-sharing model 
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where E[W] is the expected value of the wasted spaces. Here E[W] = 
6(/(p)) means that there exist positive constants c, d such that c'f(p) < 
E[W) < cf(p). Also when p — > 1 (the heavy traffic case) Coffman and 
Mitrani [5] obtained upper and lower bounds on E[W] in the form 



' ' * <E[W)<-^-(^-l 



2 v 1-p ~ j -i-pV 6 

A related model, the M/M/oo queue with ranked servers, has been 
studied by many authors Q, 0, H, [H, (12, (III. This differs from the 



current model in that if there are a total of N = Ni + N 2 spaces occupied, 
the total service rate is pN, as each customer in the restaurant is served at 
rate p. For this model various asymptotic studies appear in (TJ, (H, IITTI . In 
particular, Aldous (U showed that the mean number of the wasted spaces 



is E[W] ~ \flp log log p as p = X/p — > oo. 

A simple derivation of the exact joint distribution of finding (resp., 
N 2 ) occupied primary (resp., secondary) spaces appears in IH3| and de- 
tailed asymptotic results for this joint distribution appear in (SI, @, while 
the distribution of Max S is analyzed in 0. In B8J Rnessl showed how 
to obtain asymptotic results for the infinite server model directly from the 
basic difference equation. Since the present processor sharing model does 
not seem amenable to exact solution, we shall employ such a direct asymp- 
totic approach here. 

In this paper, we study the joint probability distribution of the numbers 
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of occupied spaces in the PS model, letting n(k,r) = Prob[Ni — k, N 2 — r] 
in the steady state. From n(k, r), we can get the probability distribution of 
the wasted spaces W from the following relations. 

oo 

Prob[w = o] = X>0'>°;i), (i.i) 

3=0 

oo 

Prob[W = L] = J2[«U,0;L + j)-7r(j,0;L + j-l)}, L > 1. (1.2) 

3=0 

Then the mean number of the wasted spaces is 

oo 

E[W] =J2 LProb l w = L \- 

L=l 

The paper is organized as follows. In section 2 we state the problem 
and obtain the basic difference equations for n(k, r). In section 3 we sum- 
marize our main results. Section 4 contains the exact solutions f or m = 1 
and m — 2, and a sketch of the derivations. In section 5 we develop a 
semi-analytic and semi-numerical method for arbitrary m, and show that 
this regains our previous results for m = 1,2. Asymptotic solutions for 
fixed m when p j 1 are obtained in section 6. Some numerical studies and 
comparisons appear in section 7. 

2 Statement of the problem 

We let Ni(t) (resp., N 2 (t)) denote the number of primary (resp., secondary) 
spaces occupied at time t. The joint steady state distribution function is 

n(k, r) = n(k, r; m) = lim Prob[Ni(t) = k, N 2 (t) — r], < k < m, r > 0. 

t^oo 
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Let p = A//i be the traffic intensity and we assume the stability condi- 
tion p < 1. The pair (iVi, iV 2 ) forms a Markov chain whose transition rates 
are sketched in Fig.l. The state space is the lattice strip {(k, r) : < k < 
m, r > 0} and the balance equations are 

(li[fc+r>o] +P) 7r(A;,r) 

+ 1 

= prt(k- 1, r)/rjfc>i] + — ■ — ir(k + 1, r)/rjfe <m ] 

A; + r + 1 

r + 1 

+ i— — 7 n(k, r + 1) + p ir(m, r - l)h k =m, r>i] • ( 2 -l) 

Here I is an indicator function. The normalization condition is 

oo m 

J^7r(M = l- (2-2) 

r=0 fc=0 

From our viewpoint we will need to consider explicitly the boundary 
conditions inherent in (|2.1[) , so we rewrite the main equation as 



k + 1 r + 1 
(1+p) 7r(fc,r) = p7r(fc-l,r) + — — 7r(k + l,r) + — — vr(/c,r + l), 

0</c<m, r>0, A; + r>0. (2.3) 



and the boundary condition at k = m is 

r + 1 

irOm, r) — p ir(m—l,r)-\ ir(m, r + 1) +p 7r(m, r — 1), r > 1. 

m + r + 1 

(2.4) 

There are also the two corner conditions 

P 7r(0,0) = tt(1, 0) + tt(0, 1) (2.5) 
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and 



(1 + p)vr(m, 0) = pn(m -1,0) + : -ir(m, 1). (2.6) 

m + 1 



In (|2.3|) when = we interpret 7r(— 1, r) as 0. The boundary condition at 
k = m in (|2.4|) can be replaced by the artificial boundary condition 

m + 1 

7r(m + 1, r) = p ir(m, r — 1). (2.7) 

m + r + 1 

This is obtained by extending (|2.3|) to hold also at k = m and comparing 



this to (|Z4|) . 



We note that the total number Ni + N 2 behaves as the number of cus- 
tomers in the M/M/l — PS queue, which is well known to follow a geo- 
metric distribution. Thus we have 

£ n(k,r) = (l-p)p N , N>0 (2.8) 

k+r=N 

and we can rewrite this as 

N 

^7r(iV-r,r) = (1 - p)p N , < N < m, (2.9) 

r=0 

N m 

Tr(N-r,r) = £ n(k, N - k) = (1 - p)p N , N > m. (2.10) 

r=N-m k=0 

These identities will provide a useful check on the calculations that follow. 
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Fig. 1 A sketch of the transition rates. 

3 Summary 

Since the analysis shall become quite involved and technical, we collect 
here some of the main results. 

For m = 1 the exact solution is given by tt(0, 0) = 1 — p and 
7r(0,r) = (l-p)(l+p) r 



i+p / | 
' 1- - ] du, r > 0, 



(3.1) 
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and then tt(1, r) = (1 — p)p r+1 — 7r(0, r + 1). For m = 2 we find that 



77(0,7-) = (l-p)(l + p) r -y +2 r(r + l) 
w 2 (l -w) 7- -^! -t) r " a * 



1 fl 



./0 



exp 



:i+p) 2 



Tt(l-U) 



X 



cfu (it, 



[l + p(l- M )(l-t)] r + 2 
7r(l,r) = (l-p)(l + P rV +2 (r + l) 2 



(3.2) 



1 rl 



^0 



exp 



L(i + p) 



;t(l-u) 



Here a* = p/(l + p) 2 < 1/4, and vr(2,r) maybe computed from 7r(2,r) = 
(1 - p)p r '+ 2 - tt(1, r + 1) - tt(0, r + 2). Also, tt(0, 0) = 1 - p, tt(1, 0) is given 
above (|430) and then tt(0, 1) = (1 - p)p - tt(1, 0). 

We next consider the heavy traffic limit where e = 1 — p — > + , with 
m = 0(1). Then setting r = Y/e, for F > and < k < m we obtain 



771 ' 

7r(fc,r) = e m - fc+1 — e - y y fc -'" 



, y m 2 + (2 — fclm — 
I - -.-( - + m + 



+ 0{e' 



This is a two-term asymptotic approximation, but it becomes invalid as 
y — > 0. In section 6 we discuss the scale y = O(e) (thus r = 0(1)) where a 
different approximation to 7r(fc, r) must be constructed. 

For moderate traffic intensities p < 1 and m = 0(1), we have 

7 777 ' 

ir(k, r) ~ (1 - p)p m +V fc - m — , r -> oo. (3.5) 

re! 

Thus the probability distribution is exponentially small (0(p m+r )), with an 
additional algebraic decay due to the factor r k ~ m . 
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4 Exact solutions for m = 



1 and m = 2 



We consider m = 1 and m = 2 for < p < 1 fixed. To solve these problems 
we use the generating function F(x, y) defined by 

m oo 

F(x, y ) = J2J2 n ^ r > k y r ( 4 - x ) 

= Zo(y)+x/ 1 (y) + ----+x m / m ( 2/ ) (4.2) 

so that 

oo 

Hy) = Y^y r < k ^)- ( 4 - 3 ) 

r=0 

Multiplying both sides of (|2.1[) by (fc + r + l)x k y r and summing over < 
fc < m, 0<r<oo we get a PDE for y): 

(p + 1 - 2px)F + (px - 1)(1 - x)F x + [(p + l)y - pxy - 

= p(y - x)x m [(m + 2)f m (y) + yf' m {y)] + 1 - p. (4.4) 

By setting x = y we find that F(y, y) — (1 — p)/(l — py). 

From (|4.1[) and (|4.2|) we have F{x,y) = fo(y) + xfi(y) for m = 1 and 



F(x, y) = fo(y) +xfi{y) +x 2 f2(y) for m = 2. Using these forms in (|4.4|) and 
comparing coefficients of x°, x 1 , and (if m = 2) x 2 , we obtain the ODEs 



m — 1 : 



h = (l + p)/ + [(i + p)j,-l]^-i + p> (4.5) 
-2pf Q -pyf = {3py - 2(1 +p)]f 1 + {py 2 - (I + p)y + l}f[. (4.6) 
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m = 2 



h = (l + p)/o + [(l+p)y-l]/o-l + P, (4.7) 
/> = (l + p)/i+ (1 + P 2 )y ~ 1 /{-p/o-^, (4.8) 

-Wi-pyfi = [^ P y-m + p)]f2 + [py 2 -(i + p)y + m- (4-9) 



We must thus solve a system of two ODEs if m — 1 and three ODEs if 
m = 2. We can simplify the calculations by observing that 

m oo oo 1 

%»)=EE' r (*' r )i^ = E^ E t(*,»-) = i -=^-. (4.io) 

fc=0 r=0 N=0 k+r=N ^ 

Therefore, for m = I, f (y) +yfi(y) = (1 - p)/(l - py) and, for m = 2, 
/o(y) + yfi(v) + y 2 J'2{y) = (1 - p)/(l - P2/)- It thus suffices to compute / if 
m = 1, and f and /i if m — 2. 

4.1 m = 1 

By solving the system of equations (|4.5[) and (|4.6|) for m = 1, we obtain 

^) = [(i +^/_ rp K 1 + ^ 2 - ( 2 + ^ - i/ksK 1 + ^x 1 - w)i + 1} 

(4.H) 

and 

A(y) = -(i-p)+ ' 



(l+p)y-l 



x 



(i +P te +T ^ + log[ '?t1 (1 T )1 - 1 V < 412 > 

1 — PS/ J 
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Here we also used the analyticity of fo(y) at the point y = 1/(1 +p), as fo(y) 
must be analytic for \y\ < 1. We first invert /o(y) and compute 7r(0, r). We 
expand (|4.11|) in Taylor series about y = 1/ ( 1 + p), setting A = y — 1/(1 + p), 
to obtain 



My) 



Replacing A in (|4.13|) by y — 1/(1 + p), using a binomial expansion of A n , 



and reversing the order of summation yields 
1-p 



fo(y) 



i + p 



1=0 n=l \ / \ ' 



+ 



(4.14) 



From (|4.14|) we can easily check that 7r(0, 0) = 1 — p and then obtain 7r(0, r) 



for r > as an infinite series: 

.(0,0 = pd-rtd+^f^-D-Q (_!_ + _£_) 

Z/ — 

(4.15) 



We can derive another form for 7r(0, r) by expanding (|4.11[) around y = 
instead of y = 1/(1 + p), which yields 

oo oo oo oo 

fo(y) = (i-p){EE(i+p) A V +3 -[2+/)+iog(i+p)££(i+p)^!/ t+J+1 

h=0 j=0 h=0 j=0 

OO OO OO OO OO 7 

ft=0 j=0 h=0 j=0 1=1 
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Therefore, the coefficient of y r when r > is 



7r(0,r) = (l-p)(l + p) r - l {(r + l)(l + p)-r[2 + p + ]og(l + p)] 



r-l 



r — I ( p 



1=1 



I \l + p 



}■ 



(4.16) 



Furthermore we have 



00 1 / 

rlog(l +p) - P + p r (l +pf- r ( i 



P 



which when used in (|4.16|) gives the following expression for 7r(0, r) when 
r > 0: 



vr(0,r) 



1-P 
1+P 

1-P 
1+P 



P 



oo / \ J— r' 

r / p 



+ P 



E 

n=0 



r + n \ 1 + p 



(4.17) 



Setting (|4.15|) equal to (|4.17|) leads to the identity 



' n + r V 1 + P 

n=l x ' 



p(l + p) r yv p) x(L + r)! 



L=0 



L! 



+ 



L+r+1 L+r+2 



which can be verified directly. We note that (|4.17|) is much more useful 



than (|4.15|) for obtaining 7r(0, r) for r — > oo. 

To get 7r(l, r), we can use either (|4.12|) , or (|4.17[) and the condition 7r(0, iV) + 
tt(1, iV - 1) = (1 - p)p N , N > 0. We thus obtain 



vr(l,r) 



(l-p)p r+1 -7r(0,r + l) 

P 



1 + P ' ' 



Si 



p 



L=0 



L 



r + L + 1 \1 + p 



L ^r + L \l + p / 
12 



(4.18) 

L 



(4.19) 



If r = 0, from (l4l9l) we get tt(1, 0) = (1 - p) J2T=i L- l p L {l + p)- L = (1 - 
p) log(l + p). This is the same as /i(0) from (|4.12|) . 

We can also derive, from (|4.17|) and (|4.19|) , alternate integral represen- 
tations for 7r(0, r) and 7r(l, r). If we write l/(r + n) in (|4.17|) as J" u n+r ~ l du, 



(|4.17|) becomes 

vr(0,r) 



(1 " P)P 



r+l 



It 



o (1 + p-up) 



-du, r > 0. 



(4.20) 



In the same way, we can rewrite (|4.19|) as 



7r(l,r) = (l-p)/ +1 (r + l) 



u' 



-du, r > 0. 



(4.21) 



1 + p - up 

Other integral representations for 7r(0, r) and 7r(l,r) can be derived from 
(|4.13|) , by replacing (n + by J" £ n dt. Then /o(y) becomes 

/o(l/) = TT-^l E(l+pt)[p(p+l)Atrp(p+l)^+ r _/ ' 



(1 + P) 2 
1-p 



n=0 



1+P 

1+P 

1-P 

1 + P 



P 

1 + 



1 + P* 



o 1 - p(p + l)yt + pt 
l+f) u 



dt + 1 
du 



l + l> n (l + P) n 



M - (1 + P)( M - 1)?/ 

l+p (u-l) 



P 
1 

A = y-— - 
1 + P 

, 1 + pt = u 



n=0 



-du 



(4.22) 



From the integral form (|4.22|) we can check that vr(0, 0) = 1 — p, and then 
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for r > we obtain (|3.1|) and (|4.18|) yields 



7r(l,r) = {1- P )p r+1 -{1- P ){l+ P y J' +P (l- 1 ^ + du 

= (i-p)(i + p y [(i + p)£ +Pr -±±(i-- \ du 

l+p / l\ r+1 

1 ) du ] , 



'U 



(4.23) 



where we used 



ti- 



ll 



l + p 



After integrating by parts in the second integral in (|4.23|) we obtain 



vr(l,r) = (l-p)(l+p) r (r + l)^" ) du, r > 0. (4.24) 



We conclude the analysis of m = 1 by giving asymptotic expansions 

for ir(k, r) as r — > oo : 



7r(0,r)~(l-p) 



7r(l, r) ~ (1 — p)p r+ , r — > oo. 



These follow easily from (I4l7|) and d419> , or (EL201) and (EL2T) . or ^TTJ and 



Km. 
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4.2 m = 2 



For the m = 2 case we can use (|4.7[) and (|4.8[) to express /2 and fi in terms 
of /o, and then obtain from (|4.9|) the following second order ODE: 



(4p 2 + 7p + 4)y 2(1 +p) 



/o(l/) 



+2 



1 + 2 

+ " — + P + P + 1 

r 2/ 



/o(y) = 2(l-p) 



1 1 1 

— - + - + l + p 

r(i - py) y 



(4.25) 

After some transformations (|4.25|) can be converted into an inhomoge- 
neous hypergeometric equation. We can show that both of the solutions 
to the homogeneous problem fail to be analytic at y = 1/(1 + p), and that 
there is a unique particular solution that is analytic at this point. To obtain 
this solution it is best to simply expand fo(y) in Taylor series, by setting 

oo 

/o = i - p + y^A 

where 



„n— 3 



(4.26) 



n=3 



;i + p) 3 a _ (i + p) 3 



p 



p 



y 



i + p. 



(4.27) 



Using (|4.26|) in (|4.25|) we can compute the A n recursively to get 



A, 



2(„-l)( n -2)(l-p)(ii#£( r 

l=n ^ 



+ pj /(/-!)(/ -2) 



1 



x- 



(n - 1 - a*)(n - a*) ■•■(/- 2 - a*)(Z - 1 - a*) 
where a* = p/(l + p) 2 . If we replace z in (|4.26|) by (|4.27|) and use the 
identity 

= E(") (- 

i=0 v 7 
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i + p 



i + p 



we can write f as 

2(l-p)(l + p) 4 



My) = i-p 3 



i=l n=i—l 



a* J (i - l)\(n - i + 1)!' 



x 

Z=n+3 



oo , v 2Z j j 

S (iTW /f/-i)f/-2)r(/-a*)- (428) 



+ py /(/-i)(/-2)r(/-a*) 

Therefore, we find that 7r(0, 0) = 1 — p and, for r > 0, 

2(l-p)(l + p) 4 



7r(0,r) 



p3 



gj" +2 ) |r ^ +2 - fl ')fe) \r ( :g-:: ) i). 



n=r— 1 



Si 



i=n+3 



2/ 



p v 1 1 



+ p7 /(/- l)(/-2)T(/-a*) 



2(l-p)(l + p) 3r+2 ^(L + r + l) 



x 



Z=L+r+2 

where L = n — r + 1. 



p r+2 L!(r-1)! v 

-iV f f-^V' ! L_ ,4.29) 

a*/ .^Al + p) l(l-l)(l-2)T(l-a*) V 7 



In a similar way we obtain n(l,r). We write as a series in ?/ using 
(|4.7[) and (|4.28|) and identify the coefficient of y r , thus obtaining 

tt(1,0) = P(1-P) + 2(1 ~ P) 1 1 + P)3 [f:n(n + l)(-i) n r(n + l-a 

/ P \ 1 1 

X 

l=n+2 



Si 



+ p7 /(/- l)(/-2)T(/-a*) 
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and for r > 

2(r + l)(l-p)(l + p) r+5 v 



vr(l,r) 



oo 



E 



x r(» + 2-„')E (^) 



n=r— 1 
2i 



(n + l)(n + 2)!(-l) 
(n — r + 1)! 

1 1 

/(/-l)(/-2)r(/-a*)' 



n— r+1 / 

a* 



(4.30) 



From (|4.29|) and (|4.30|) we can obtain alternate integral representations 



for 7r(0, r) and 7r(l, r). We first rewrite (|4.29|) and (|4.30|) as 

(L + r + 1)! / p 



oo oo 



vr(0,r) = ^l-rffl+rV^E 



L=0 j=0 



L!(r-1)! \l + p 



x 



vr(l,r) 



l T(L + r + l-a*) T(j + 1) 
T(j + L + r + 2-a*) j\ 
1 

(j + L + r + 2)(j + L + r + l)(j + L + r) ' 

2 ( r + i)(i - „)(! + ,r y« f E (L + r)( ^,: r + 1)! 

2j 



(4.31) 



L=0 j=0 



p x - , L r(L + r + l-a*) r(j + l) 



1 + P 



{-pY 



T(j + L + r + 2-a*) j\ 
1 



(j + L + r + 2){j + L + r + 1)(J + L + r) 
We use the beta function 



(4.32) 



B(j + l,L + r + l-a*) 

_ r(j + l)r(L + r + l- q*) 
T(j + L + r + 2 - a*) 



L+> — u* 



tft (4.33) 
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and represent the last factors in (|4.31|) and (|4.32[) as the following integral 

1 

(j + L + r + 2){j + L + r + l)(j + L + r) 

1 1 1 

+ 



2(j + L + r + 2) j + L + r + 1 2(j + L + r) 



1 



2 Jo 



(1 - u) 2 u 1+L+r - 1 du 



u 



'(l-u) j+L+r - l du. 



(4.34) 



Using (14331 and d434l) in (14311) and (1432b leads to 



vr(0,r) = (l-p)(l + P rV +2 f fu\l-ur\l-ty- a * 

Jo Jo 



oo oo 



X 



X 



L=0 j=0 



L!(r - 1)! 



1 -{l-u) j dudt, r>0 



(435) 



and 



vr(l,r) = ( r + l)(l-p)(l + pr-y+ 2 /' f u\l - u) r -\l - t) r ~ a ' 

Jo Jo 



oo oo 



££ ' t + r "^ r+1|! N'[(i-)H-f 



L=0 j=0 



L! r! 



x 



,1 + P/ J' 
Using the binomial series 



'' » -(l-tt) j (itt^ r>0. 



(436) 



L=0 



(L + r + 1 ) ! z L = (r + 1) , (1 _ z) -,-2 



(437) 



we can explicitly evaluate both infinite series in (|435|) and (|436|) , to obtain 
the integral representations (B.2[) and (|33|) for 7r(0, r) and 7r(l,r), for r > 
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0. We can similarly obtain 7r(2, r) by using (|4.8[) , but we can also use the 
identity vr(2, r) = (1 - p)p r+2 - tt(1, r + 1) - tt(0, r + 2) which holds for all 
r > 0. Also, 7r(l, 0) can be computed from p(l — p) — 7r(0, 1). 

We again conclude by giving asymptotic formulas for 7r(/c, r) for r — > 
oo. These follow easily from the double integrals in (|3.2|) and (|3.3|) , since 
for r — > oo the integrands become concentrated in the range u, t — 0(l/r). 
Therefore, evaluating the integrals by the Laplace method we find that as 

r — > oo 

T(0,r) ~ ^(1-P)p r+2 , 

vr(l,r) ~ -(l-py +2 , 
r 

7r(2,r) ~ (l-p)p r+2 . 

5 Semi-numerical method for fixed m 

In this section we discuss how the exact solutions for fixed m can be ob- 
tained by a semi-numerical approach. We shall reduce the solution of the 
two-dimensional problem in (|2.1|) to a one-dimensional one, which must 
be solved numerically for general m. 

If we set ir(k, r) = D(N, r) where N = k + r, (|2.3|) becomes 

(N + 1) [(1 + p)D(N, r) - pD(N - 1, r) - D(N + 1, r)] 

= (r + l)D(iV + 1, r + 1) - rD(N + 1, r). (5.1) 

Since (|5.1[) is separable, we seek solutions in the product form D(N, r) = 
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a(N)(3(r). Then (|57T) becomes 

(N + l) 



, a(iV) a(JV-l) 
1 + P) 



r + l)^±l)-r = -M 

(3(r) 



a(JV+l) a(iV + l) 

where — M is a separation constant, which we take to be a negative integer. 
Solving the difference equations for a and (3 we obtain 

«(A;+r,M) = — / / - w (l-zrW-dfl-pz)^ 1 "')^, (5.2) 

2vn J c (1 - z)(l -pz) 

where M < k + r and C is a small loop around z = 0, and 

^• M >= ,. (M-r)! ' M£ r ' 
Therefore, we can write 7r(fc, r) as a sum over all r < M < k + r. 

n(k, r) = ^ C(M) J _ a(fc + r, M), r>0 (5.3) 

where C (M) is a function of M. 

The constants C(M) will be chosen to satisfy the boundary condition 
along k = m, as given by (|2.4|) . We are allowing 7r(fc, r) to have a different 
form for r = 0, and we shall compute ir(k, 0) from (|2.3|) (with r = 0), (|2.5[) , 
and fl53). We define C(0) in (O by requiring (l53|> to hold at r = if 
k = m. Thus (|5.3|) will apply for {0 < k < m, r > 1} or {k = m, r > 0}. 

When k = 0, (|5.2[) yields a(r, r) = 1 and if = 1, we obtain a(r + 1, r + 
1) = 1 and a(r + l,r) = (1 + p)(r + 1). Note that a{N, M) = if M > N. 
Thus (|5.3[) yields 

vr(0,r) = ^, r>0, (5.4) 
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and 



tt(1, r) = + (i + p )(i + r )£M, r > 0. 



7*! T\ 

If we set C(r) = r\d(r), (|5.3|) can be written as 



(5.5) 



7r(fc, r) 



Y^(-l) M ' r d{M)( )a{k + r,M) 



M=r 
k 



J2(-l) l d(r + l)( l + r )A(L,r:l). 



(5.6) 
(5.7) 



for r > or k = m, where M = r + I and 



A(k,r;l) 



1 



J-k-l 



2m J c (1 — z)(l — pz) 
From (|5.8[) we obtain 



(1 _ z )-P+r)/(l-p) (1 _ p^+O/d-p)^. 

(5.8) 



A(Z,r;Z) = 1, 
A(Jfe + l,r-l;Z + l) = A(Jfe,r;Z). 



(5.9) 
(5.10) 



We solve for n(k, 0) in terms of ir(k, 1). If we define A(fc) =p fc [7r(A;,0) — 
7r(A; — 1, 0)], (|2.3|) with r = can be written as 

tt(M) 



A(Jfe) - A(Jfe + 1) 



(k + l)p fc + 1- 
Summing (|5.11|) over 1 < /c < k — 1 gives 

From the corner condition (|2.5[) and (|5.12|) we get 



(5.11) 



(5.12) 



fc-i 



7T 



(*, 0) - n(k - 1, 0) = -p k ~\l - P f-J2 l> J 



j=0 
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Summing over 1 < k < k we obtain 



fc-l 3 



7i(k, 0) 



3=0 1=0 
fc-1 



Z + l 

1 - p k - 



1=0 



I 1-P J 



(5.13) 



This holds also if k = 0, since then the sum is void and we obtain vr(0, 0) = 
l-p. 

We can compute 7r(m, 0) (the probability that all m primary and no 
secondary servers are occupied) from either (|5.7[) or (|5.13|) (with (|5.7|) used 
to compute 1) in the sum). By equating these two expressions we find 
that 

:i - P )p m 



d(0) 



A(m,0;0) : 



(5.14) 



and also obtain the identity 



™ n (l _ m-l+l\ 

A{m, 0;n) = Y^ Th T\ 0; n) 



l=n 



1(1 - P ) 



which can also be derived directly from (|5.8[) . 

We still need to determine d(r) (and hence C(r)) for r > 1. To this end 
we apply (|2.4|) (or (|2.7[> ), which yields 

m + 1 



m + r + 



J2(-l) l d(r + l)( + )A(m + l,r;l) 
i=o V r / 

= p ^(-l)'d(r + /-l)l r JA(m,r-l;Z), r > 1 (5.15) 

and then we must also satisfy the coiner condition in (|2.6|) . We can view 
(|5.15|) as a difference equation for the d[r), of order m + 2. In general this 
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must be solved numerically, though below we discuss the cases m = 1, 2, 3. 
For general m we have thus reduced the two-dimensional problem (|2.3|) to 
the one-dimensional one in (|5.15|) . 



5.1 m = 1 



When 771=1, from (|5.15|) we obtain the following third order difference 
equation: 

2(1 + p)(r + 1) d(r + 1) - (r + 1) d(r + 2) + p(l + p)r d(r - 1) 

+ [p(l -r) - (1 +p) 2 (r + 1)] d(r) = 0, r > 0. (5.16) 



Rearranging (|5.16|) gives 



p(l + p) [(r + l)d(r) - rrf(r - 1)] - (1 + 2p) [(r + l)d(r + 1) - rd(r)] 

+ (r + l)d(r + 2) - rd{r + 1) - [d(r + 1) - d{r)\ = 0, r > 0. (5.17) 

If we define i?(r) = p(l + p)(r + l)d(r) - (1 + 2p)(r + l)d(r + 1) + (r + 
l)d(r + 2) - d(r + 1), then (f5l7) yields 7?(r) = i9(r - 1), so that 0(r) = i?(0) 
is a constant. To find this constant we use the corner condition (|2.6|) . From 
either (|5.5|) (which holds now also at r = 0) or (|5.7[) and (|5.8|) , (|2.6|) can be 
rewritten as 



p(l - p) = (1 + p) 2 rf(0) - 2(1 + p)d{\) + d{2). 



(5.18) 
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Since d(0) = p(l - p)/(l + p) by d5T4l) , (l5l8l) implies that $(0) = p(l + 
p)d(0) - 2(1 + p)d(l) + d(2) = 0. Thus #(r) = for all r > and hence 

0(r) = p(l+p)(r + l)d(r)-[(l + 2p)(r + l) + l]d(r + l) + (r + l)d(r + 2) = 0. 

(5.19) 

We solve (|5.19|) by using the generating function 



r=0 



Multiplying (|5.19|) by z r and summing over all r we obtain 

z[l - (1 + p)z]G' - [1 + (1 + p)z\G = -d(0) : 1 : 



1 — pz 

The solution to (l5T20> is, using d(0) = p(l - p)/(l + p) = G(0), 
P(l-P) l\ 1+P 



(5.20) 



G(z) 



(l + p)[l-(l + p)zP 



z\og[(l-pz)(l + p)]-(l + p)z 



(5.21) 



The integration constant that arises in solving (|5.20[) was determined by 
using the analyticity of G(z) at z = 1/(1 + p). 



We compare G{z) = Y,Z Q d{r)z r to /o(*) = E^cMM^ in SI). 
Since, for r > 0, 7r(0,r) = d(r) from (|5.4|) , the coefficients of G(z) and 
/o(^) differ only when r = 0. Then we use 

G(0) = d(0) = p^4 = /o(0) P 



;i + p) — i + p 

to find that f (z) - G{z) = (1 - p)/(l + p) = vr(0, 0) - pn(0, 0)/(l + p) 
tt(0,0) -d(0). 
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5.2 m = 2 



The boundary condition (|2.4[) for /c = m = 2 is 



(1 + p) tt(2, r) = p tt(1, r) + tt(2, r + 1) + p tt(2, r - 1). (5.22) 



r + 3 



From (|5.7|> - (|5.10|) we obtain 



vr(l,r) = d(r)A(l,r;0)-d(r + l)(r + l) 

= d( r )(l + p)(r + l)-d(r + l)(r + l), (5.23) 



and 

vr(2,r) 



d(r)A(2, r; 0) - d{r + l)(r + 1)A(2, r; 1) + d(r + 2) 



+ l)(r + 2) 



d(r)[(l + p + p 2 ) + (1 + p)V]^^ 



d(r + + p)(r + l)(r + 2) + d(r + 2) 



(r + l)(r + 2) 



(5.24) 
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Using (|5.23|) and (|5.24|) in (|5.22[) yields the following fourth order difference 
equation 

[d{r + 3)(r + l)(r + 2) - d(r + 2)r(r + 1)] 

- 3(1 + p)[d(r + 2)(r + l)(r + 2) - d(r + l)r(r + 1)] 
+ [d(r + 2)r(r + 1) - d(r + l)r(r - 1)] 
- p(l + p + p 2 )[d(r)(r + 2) - d(r - l)(r + 1)] 
- p(l + p) 2 [d(r)r(r + 2) - d(r - l)(r - l)(r + 1)] 

+ (1 + p + p 2 )[d(r + l)(r + 3) - d(r)(r + 2)] 
+ (1 + p) 2 [d(r + l)(r + l)(r + 3) - d(r)r(r + 2)] 
+ p(l + p)[d(r + l)r(r + 2) - d(r)(r - l)(r + 1)] = 0. (5.25) 

If we define 

tf(r) = d(r + 3)(r+ l)(r + 2) - 3(1 + p)d(r + 2) (r + l)(r + 2) 

+ d(r + 2)r(r + 1) - p(l + p + p 2 )d(r)(r + 2) - p(l + p) 2 d(r)r(r + 2) 

+ (1 + p + p 2 )d(r + l)(r + 3) + (1 + p) 2 rf(r + l)(r + l)(r + 3) 

+ p(l + p)rf(r + l)r(r + 2), (5.26) 



then (|5.25|) becomes i?(r) = ?9(r — 1), so that -#(r) = $(0) is a constant. From 



(|5.14|) we obtain d(0) for m = 2 as 



= 7~T~~T (5-27) 

1 + P + 
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To find the constant $(0) we again use the corner condition (|2.6|) . Using 
(EH, d527b , (15T31) , and (JHTJ, (O for m = 2 becomes 



d(3) = p 3 (l - p) - -(6 + 9p + 6p 2 )d(l) + 3(1 + p)d(2). (5.28) 



Using (15T281) in d5T26l) we find that i?(0) = and thus tf(r) = for r > 0. 

Now the coefficient functions in (|5.26|) are quadratic in r, and we intro- 
duce again the generating function G(z) = YlT=o d(r)z r to solve i?(r) = 0. 
This leads to the second order ODE 



[l-2z- 3pz + (1 + p)(l + 3p)^ 2 - p(l + p) V ] G 



+ 



2 + (6p 2 + 9p + A)z - p^ 2 (4p 2 + 7p + 4) 



G' 



+ 



4 + - + 2-2p(p 2 + p + l)^ 



2r 



G 



P 2 (l-P) 



1 + p + p 2 V £ 2 



1 1 
+ - + 1 

z 



(5.29) 



Since G(z) = f (z) - tt(0, 0) + d(0) = f (z) - (1 - p) + p 2 (l - p)/(l + p + p 2 ), 
if we write (|5.29|) in terms of / 0/ we obtain (|4.25|) (after dividing (|5.29|) by 
1 — pz). Equation (|5.29|) can be solved as in section 4, to ultimately obtain 
d(r), and then %(k, r) from (|5.7|) . 



5.3 m > 3 



This method can be used for general m, but we find that it becomes very 
difficult to solve the ODE corresponding to (|5.20[) or (|5.29|) when m > 2. 
For example, when m = 3 we obtain from (|2.4[) or (|5.15|) a fifth order dif- 
ference equation for d(r). This equation is a perfect difference that may be 
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summed to yield the fourth order difference equation 

(r + l)(r + 2)(r + 3)d(r + 4) + [r(r + l)(r + 2) 

- 4(1 + p)(r + l)(r + 2)(r + 3)]d(r + 3) + 3[(2 + p)(r + l)(r + 2) 
+ (1 + p)(l + 2p)(r + l)(r + 2)(r + 3) + 2(r + l)]d(r + 2) 
+ [3(2 + p)(r + l)(r + 2) - 12(1 + p)(l + p + p 2 )(r + l)(r + 2) 
- (1 + p) 2 (l + 4p)r(r + l)(r + 2) - 2(3 + 2p)(2 - p)(r + 1) + Qr}d(r + 1) 
+ (1 + p)[p(l + p) 2 (r - l)r(r + 1) + 9p(l + p + p 2 )r(r + 1) 

- 3p(l + p 2 )(r - 2)(r + 1) + 2p 2 r]d(r) = 0. (5.30) 

Multiplying (|5.30[) by z r and summing over r we obtain the following third 
order ODE for G(z) for m = 3: 

[ p(l + p) 3 z 3 - (1 + p) 2 (4p + l)^ 2 + 3(1 + p)(2p + l)z - 4(1 + p) 
+ M + ^ ] G'" + 3 [ p(l + p)(3p 2 + hp + 3)z 2 

-(8p 3 + 17p 2 + 13p + 3)z + 2(3p 2 + hp + 3) - (l + i + ] G" 
+2 [ p(l + p)(9p 2 + lOp + 9)z - 2(1 + p)(6p 2 + hp + 6) 
+3 fl + i + i + ] G' + 6 [ p(l + p)(l + p 2 ) 

f 1 , 1 , 1 , 1 W_ 6p 3 (l-p) A J_ 1 JL\ 

% 2 2 3 (i + p)(i + p 2 ) ^ 2 ^ 3 ' 

(5.31) 

This appears too difficult to solve in terms of elementary or special func- 
tions. Thus, for general m we shall obtain n(k, r) only asymptotically, as- 
suming that p t 1- 
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6 Asymptotic solutions for p f 1 with fixed m 

In this section, we fix m and let p — > 1. We thus set p = 1 — e where e 
is small and positive. We shall analyze the problem for two ranges of r, 
r = Y/e = 0(e~ l ) and r = 0(1). Note that necessarily k < m = 0(1). 
In this heavy traffic limit most of the probability mass occurs in the range 

r = 0(s- 1 ). 

6.1 r = 0(£ _1 ) 

First we analyze (|2.3|) for r = 0(e _1 ) by setting 

r = e- 1 Y, 7r{k,r) = e m ' k n{k,Y;e). 
On this scale (|2.3|) becomes 



(2 - e)[Y + e(k + l)]TZ(k,Y) = e(l — e)[Y + e(k + l)]1Z(k — 1, Y) 
+(k + l)TZ(k + 1,Y) + (Y + e)TZ(k,Y + e), 

while the boundary condition (|2.4|) at k = m becomes 

(2 - e)[Y + e(m + l)]K(m, Y) = (Y + e)7Z(m, Y + e) 

+ (l-e)[Y + e(m + l)][K(m,Y -£)+eK(m-l,Y)\. (6.1) 

We expand TZ(k, Y) = TZ(k, Y; e) in the form 

TZ(k, Y) = e[TZ i0) (k, Y) + eK {1 \k, Y) + 0(e 2 )}. 
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From (|2.3[) and (|2.7|) we obtain the following equations for IZ^ and TZ^\ 
for y > 0, 

Y 1l i0) (k,Y) = (k + l)n i0 \k + l,Y), 0<k<m (6.2) 

y ^-n (0) (k,Y) + (y - 2k - 1) n {0) (k, y) + y^ (0) (k -i,y) 
yy 

= y^ (1) (A;,y)-(A; + l)^ (1) (A; + l,y), < k < m. 

(6.3) 

The general solution to (|6.2|) is 

^(°) (Jfe y) = ^-^ (0) (0, y ) , < fc < m, Y > 0. (6.4) 
k\ 

Using (|6.4|) in (|6.3|) and solving for T^A 1 ) yields 



Y k „ n , y 



^W(A:,y) = — ^ (1) (o,y) 



fc-i 



fc! v ' ' {k-l)\ 



(y - i)^ (0) (o, Y) + y A^(o)( 0; y) 



(6.5) 

It remains to determine (0, Y) and (0, Y) . We use the boundary 
condition at k = m. From (|6.1|) we find that this equation is satisfied au- 
tomatically to orders 0(1) and 0(e), while at orders 0(e 2 ) and 0(e z ) we 
obtain 

' r 7^ (m, y ) + t^7^ (0) (m, y ) = 0, y > (6.6) 



<9y 2 v ' ' dY 
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and 

Y ^n (1 \m, Y) + (Y- m)^-n {1) (m, Y) + F^ (1) (m - l,Y) + ft (1) (m - 1, F) 

__J_^(o) ( m _iy)_ F^-^ (0) (m - 2, Y) + n (0) (m — 1, F) — ft {0) (m — 2, Y). 
oY oY 



Since ^ (0) (m, Y) = ^ (0) (0, Y) Y m /m\ should be integrable about Y = oo, 
from (|6.6|) we conclude that 7^°) (m, F) is proportional to e~ Y , so that 



and hence 



rt°\0,Y) = ^e- Y h(m), (6.8) 



7r(fc,r) s rn - k Y k - m e- Y h(m) 
k\ 



for some function h{m). The normalization condition (|2.2|) then gives h(m) 
1. From (|6.4[) and (|6.7[) we obtain the following equation for "T^ 1 ) (0, Y) 



Y ^r^\0,Y) + (F + 2m)A^ 1 )(0,F)+(^-^ + m)^ 1 )(0 ) F) 

!)w<°>(o,y) 



2m 2 + 4m m 2 + 2m F " 

F 2 + F~ 
2m 2 + 4m m 2 + 2m F\ m! _ Y 



y2 + — ijy^ <6 ' 9) 

Solving the ODE (|6.9|) we find that the general solution takes the form 

F 2 \ m\ 



K {1) (0, F) = Ci(m)F- m + C 2 (m)Y- m e- Y - (m 2 + 2m + 



(6.10) 



(6.7) 
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To fully determine H^(0,Y) we use (|2.10[) , which can be written as, for 

N >m, 



e(l-e) N = ^ n(m — L, N + L — m) 

L=0 

m 

= ^ E L H(m - L,eN + eL - em) 

L=0 

Setting N = Z/e and using 

ry "1 

(l- £ ) z/£ = exp -log(l-e) =e- z 1--Z + 
& - 2 



(6.11) 



from (|6.11|) we find that 



n {0 \m,z) + e 



d 



n {1 \m, Z) - m—n {0 \m, Z) + 1l {0) (m - 1, Z) 
uZ 



e- z [l--Z + 0(e% 



(6.12) 



By comparing terms for order 0(1) in (|6.12[) we obtain lZ^(m, Z) = e z , 
which is the same as what we obtained from (|6.4|) and (|6.8|) . Using (|6.4|) , 
(1631) , d£5b and (f6TT2b , the 0(e) terms in d6Tl2b yield 



n {1 \o,z) 



-e e 



Z m 2 + 2m 

h m H 

2 Z 



(6.13) 



Therefore, by comparing (|6.10[) and (|6.13|) we find that C\ — and C2 = 
-m m!in (|6.10|) , and thus obtain (|3.4|) , which is a two- term asymptotic ap- 
proximation to n(k, r) for r = 0(e _1 ). Note that the second term becomes 
comparable to the leading term for small Y, with Y = 0(e). 
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6.2 r = 0(1) 

We next examine the problem for r = 0(1). Now both k and r are 0(1) 
and we set ir(k, r) = eQ(k, r; e), which satisfies 

(2-e)(k + r + l)Q(k,r) = (k + l)Q(k + 1, r) + (r + l)Q(k, r + 1) 

+ (l-£)(k + r + l)Q(k-l,r) (6.14) 

for < k < m, r > 0, with the boundary condition 

(2-e)Q(m,r) = (1 - e)Q(m - 1, r) + (1 - e)Q(m, r - 1) 

+ + r>l, (6.15) 

m + r + 1 

and the corner conditions 

(2-e)(m + l)Q(m,0) = (1 -e)(m + l)Q(m - 1,0) + Q(m, 1), 

and 

(1 - e)Q(0, 0) = Q(l, 0) + Q(0, 1). (6.16) 
We expand Q(k, r) in the form 

Q(k, r) = Q<-°\k, r) + eQ {1 \k, r) + 0(e 2 ). 



Then from (|6.14|) we obtain at the first two orders 



2 Q^\k,r)-Q^\k-l,r) = _A±1- gW(A ; +l, r ) + ^±L. QW^.r+l), 

/e + r + 1 /e + r + 1 

(6.17) 



33 



and 



2Q (1) (fc,r) -Q (1 \k- l,r) 



1 h 1 Q«(]fe + 1, r) - , r + 1 - Q«(]fe, r + 1) 



k + r + 1 
Q ( -°\k,r)-Q ( -°\k-l,r 



k + r + 1 



(6.18) 



Equations (|6.17|) and (|6.18|> hold for < < m, r > 0, + r > 0, with 
Q(— 1, r) = 0. We also obtain boundary conditions at k = m from (|6.15|) as 

' 1 Q(°)(m,r+1) (6.19) 



(m,r)—Q { } (m — l } r) — Q ( >(m } r — 1) 



)(0) 



(0), 



and 



m + r + 1 



r + 1 



2Q w (m,r) - Q ll) (m - 1, r) - Q U) (m, r - 1) - 

m + r + 1 

= Q (0) (m, r) - Q (0) (m - 1, r) - Q (0) (m, r - 1) 



g (1) (m,r + 1) 



The corner conditions at (k, r) = (0, 0) are given by, in view of (|6.16|) , 



Q(°)(1,0) + Q(°)(0,1) = 1 

q«(i,o) + q«(o,i) = -l 



(6.20) 
(6.21) 



where we used Q (0) (0, 0) = 1 and Q (1) (0, 0) = 0. 

We also require that the expansions on the Y and r scales asymptoti- 
cally match, in an intermediate limit where Y — > and r — > oo. This means 
that 

e m " fc [ ft (0) (A;, F) + (Jfe, y ) + 0(£ 2 ) ] |y^ 

~ Q (0) (A:, r) + eQW (k, r) + 0(e 2 ) I™ . (6.22) 
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By setting Y = er in the left side of (|6.22[) and using (|3.4|) we find that, for 

r — > oo, 



Q {0) (k,r) 



mi 
~k\ 



„k—m 



\ _ m 2 + {2- k)m - k + Q , _ 2 " 



(6.23) 



and 



QW^, r) = ^ r fc ~ m [-r + (m - k)(m + 1) + 0(r -1 )]. (6.24) 

We also get E k+r =N Q (0) (k,r) = 1 and £ fc+r=JV Q (1) (M = -iV from 

The basic problem (|6.17|) , (|6.19|) , and (|6.20[) for the leading term (A;, r) 
is only slightly simpler than the full problem in (|6.14[) (or (|2.3|) ). The only 
simplification is that p has been replaced by 1. We have not been able to 
solve this for general m, but the results in section 4 can be used to identify 
Q^(k,r) and Q ^'(k, r) when m = 1 and m = 2. For m = 1 we use (|3.1|) 
and (|4.23|) and expand these expressions for p — > 1 and r > 0, to find that 



Q (0) (0,r) 
Q«(0,r) 
Q(°)(l,r) 

gW(l,r) 



2' 



r-l 



-(r - 1)2 

»2 



1 



r-2 



li 



1 



du — 



1 - 2 r I (1-^1 ^ 

1-Q(°)(0,r + 1), 

-r- 1-Q (1) (0,r + 1). 



-r - - 4 
2 



r+l 



(6.25) 
(6.26) 
(6.27) 

(6.28) 
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For m = 2 (|3.2|) and (|3.3|) give the leading terms as, for r > 0, 



Q lUj (0,r) = 2 r " 2 r(r+l) / / exp 

'o Jo 



1 /•! 



g (0) (i, 



X 



2 r (r + lj" ; / exp 
'o Jo 

r-2(l -u)(l-t) 



1 /•! 



t(l-u)' 

4 



u 2 (i - u) r -\\ - ty-v* 



(2-t-u + ut) 



r+2 



du dt 



(6.29) 



u 2 {\ -uy-\\ -ty- 1/4 



(2-t-u + uty+ 3 

and (2, r) can be computed from 



du dt. 



(6.30) 



Q(°> (2, r) = 1 - Q(°) (1, r + 1) - g (0) (0, r + 2). 



(6.31) 



The results of section 4 can also be used to identify the correction term 
QW{k,r) for m = 2 and fc = 0, 1,2. From ((625T) - (|6311) we can easily verify 
that (|6.23|) and (|6.24|) are satisfied for m = 1 and m = 2, by expanding the 
integrals in (|6.25|) - (|6.30|) for r — > oo by the Laplace method. 



7 Numerical Studies 



We assess the accuracy of some of the asymptotic formulas we obtained. 
First we consider the heavy traffic case, where (|3.4|) applies. In Table 1 
we consider e — 1 — p — 0.1,0.05,0.02 and 0.01. We take m = 3 and 



Y — er — 1, and compare the one and two term approximations in (|3.4[) 
to the exact (numerical) values of 7r(fe, r). For each value of e we give 
ir(k,r) for < < 3. The two-term approximation includes the 0(e) 
term inside the brackets in (|3.4[) . When e = 0.1 we see that the one-term 
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approximation is quite poor, while two-term approximation is even worse, 
and may leads to a negative answer. However, as e decreases to 0.01 the 
agreement becomes quite good, and we also clearly see the improvement 



obtained by using the second term in (|3.4|) . 

In Table 2 we consider the limit r — ► oo with k,m = 0(1), where we 
have the asymptotic formula it(k,r) ~ (1 — p)p m+r r k ~ m m\/k\. We begin 
with r = 5, which we ultimately increase to r = 50. The agreement be- 
tween exact and asymptotic results is quite poor when r = 5 but becomes 
much better for r = 50. For each value of r the agreement is the worst 
when k — and improves with increasing k. These data suggest that it 



may be useful to compute also the correction term, as we did in (|3.4|) . 

These comparisons show that the asymptotics agree reasonably well 
with the exact numerical values of ir(k, r). In some cases it proves useful 
to compute more than one term in the asymptotic series. 
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Table 1 



p — 1 — e, m = 3, Y — 1, r — Y/e. 



€ 


fc 


exact 


one-term 


two-term 


0.1 





5.40 x 10~ 5 


.000220 


< 




1 


.000696 


.00220 


< 




2 


.00463 


.0110 


< 




3 


.0210 


.0367 


.0128 


0.05 





6.28 x 10~ 6 


1.37 x 10~ 5 


1.03 x 10~ 6 




1 


.000146 


.000275 


7.58 x 10" 5 




2 


.00172 


.00275 


.00131 




3 


.0136 


.0183 


.0124 


0.02 





2.50 x 10~ 7 


3.53 x 10~ 7 


2.22 x 10~ 7 




1 


1.34 x 10~ 5 


1.76 x 10~ 5 


1.25 x 10" 5 




2 


.000361 


.000441 


.000348 




3 


.00649 


.00735 


.00640 


0.01 





1.84 x 10~ 8 


2.20 x 10~ 8 


1.79 x 10~ 8 




1 


1.91 x 10- 6 


2.20 x 10~ 6 


1.88 x 10~ 6 




2 


9.96 x 10- 5 


.000110 


9.87 x 10~ 5 




3 


.00345 


.00367 


.00343 
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Table 2 



p = 0.5, m = 3, 5 < r < 50. 



r 


k 


exact 


asymptotic 


5 





2.29 x 10~ 5 


9.38 x 10~ 5 




1 


1.58 x 10~ 4 


4.69 x 10" 4 




2 


6.02 x 10~ 4 


1.17 x lO" 3 




3 


1.65 x 10~ 3 


1.95 x 10~ 3 


10 





1.60 x lO" 7 


3.66 x 10~ 7 




1 


1.94 x 10~ 6 


3.66 x lO" 6 




2 


1.23 x 10~ 5 


1.83 x 10~ 5 




3 


5.49 x 10~ 5 


6.10 x 10~ 5 


20 





2.83 x 10- 11 


4.47 x 10- 11 




1 


6.29 x 10~ 10 


8.94 x 10~ 10 




2 


7.18 x 10~ 9 


8.94 x 10~ 9 




3 


5.60 x 10~ 8 


5.96 x 10~ 8 


30 





9.41 x 10~ 15 


1.29 x 10~ 14 




1 


3.04 x 10~ 13 


3.88 x 10- 13 




2 


5.00 x 10~ 12 


5.82 x 10~ 12 




3 


5.57 x 10- 11 


5.82 x 10~ n 


40 





4.18 x 10~ 18 


5.33 x 10~ 18 




1 


1.77 x 10~ 16 


2.13 x 10~ 16 




2 


3.79 x 10- 15 


4.26 x 10~ 15 




3 


5.49 x 10~ 14 


5.68 x 10~ 14 


50 





2.19 x 10~ 21 


2.66 x 10~ 21 




1 


1.15 x 10~ 19 


1.33 x 10~ 19 




2 


3.03 x 10~ 18 


3.31 x 10~ 18 




3 


5.40 x 10~ 17 


5.55 x 10~ 17 
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